A Fluctuation-Dissipation Process without Time Scale 
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We study the influence of a dissipation process on diffu- 
sion dynamics triggered by fluctuations with long-range cor- 
relations. We make the assumption that the perturbation 
process involved is of the same kind as those recently studied 
numerically and theoretically, with a good agreement between 
theory and numerical treatment. As a result of this assump- 
tion the equilibrium distribution departs from the ordinary 
canonical distribution. The distribution tails are truncated, 
the distribution border is signalled by sharp peaks and, in the 
weak dissipation limit, the central distribution body becomes 
identical to a truncated Levy distribution. 
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I. INTRODUCTION 

The derivation of thermodynamics from dynamics is 
still an open field of investigation 0J|] . Hereby we focus 
on a related but seemingly less ambitious purpose, the 
derivation of fluctuation-dissipation processes from de- 
terministic dynamics It has been recently pointed 
out [Q that a genuinely dynamic derivation of Brown- 
ian motion would be essentially equivalent to a mechani- 
cal foundation of thermodynamics, thereby implying that 
also this avenue might be fraught by strong conceptual 
difficulties. It has been remarked Q that the dynamic 
foundation of Brownian motion, as described by an ordi- 
nary Fokkcr-Planck equation, imply fluctuations with a 
finite correlation time t, namely, rests on the existence of 
a finite microscopic time, or, equivalently, on the micro- 
scopic foundation of the linear response theory [|| . How- 
ever, the resulting transport equation can be identified 
with a bona fide Fokker-Planck equation only if |6| the 
corresponding relaxation process is exactly, not approx- 
imately, exponential: a property in harsh conflict with 
both quantum M and classical |g[ dynamics. This is the 
main reason why the problem of the dynamic foundation 
of the ordinary Fokker-Planck equation is not yet set- 
tled and further efforts must be made not excluding the 
possibility of either non-Newtonian effects j| or sponta- 
neous fluctuations || , both implying a kind of generaliza- 
tion of ordinary classical and quantum mechanics. Here 
we reverse the perspective and rather than imposing the 
Markovian approximation, incompatible with the deter- 
ministic nature of the system under study, we discuss the 
consequence of explicitly rejecting the requirement of a 



finite microscopic time scale. To conduct this discussion, 
we adopt Occam's principle, namely we study the sim- 
plest dynamical system with the essential features neces- 
sary to produce a fluctuation-dissipation process without 
using the requirement of a finite microscopic time scale. 
Let us consider the Liouville-like equation 

d 

-^_Pt(x, t w, t) = Lp T (x, £, w, t) 

= [-^+f(-A 2 x)]p T (x,^w,t). (1) 

This means that we study for simplicity a one- 
dimensional case. The one-dimensional variable of inter- 
est x undergoes the influence of a "fluctuation", called 
£. The dynamics of £ is driven by the operator P which 
expresses concisely the action that a set of variables w 
can exert on £ so as to render disordered its time evo- 
lution. Thus, in principle, the "stochastic" dynamics of 
x can be either provoked by nonlinearity or by the large 
number of degrees of freedom. The unperturbed fluctu- 
ations £ are the source of diffusion of the variable x. To 
undergo also dissipation, the second key ingredient of a 
fluctuation-dissipation process, the variable x must also 
exert a feedback on the dynamics of £. This important 
property is expressed by the dependence of T on A 2 , left 
unspecified. As we shall see, we shall assume a linear 
departure from the unperturbed condition Tq given by 
f (-A 2 x) = f o - A 2 xf i. The operator fi drives the 
bath response to an external perturbation M. The pa- 
rameter —A 2 x denotes the strength of the feedback and 
the minus sign alludes to the reaction nature of the effect. 
Pursuing our program inspired to Occam's criterion, we 
are forced to assume the variable £ to be dichotomous. 
It would be surprising if £, and so the microscopic statis- 
tics, would be Gaussian. In a sense, there would be no 
problem to solve at all. Therefore we must adopt a non 
Gaussian statistics. Thus, we fix the statistics of £ to be 
dichotomous, since dichotomous statistics seem to be the 
simplest example of non-Gaussian statistics. 

We plan to prove that a bath, described by the un- 
perturbed operator To with a diverging correlation time, 
r = oo, yields a form of equilibrium strongly departing 
from the ordinary canonical prescription. The proof is or- 
ganized as follows. In Section II we derive a generalized 
master equation. We make the basic assumption that 
the memory kernel of this master equation only depends 
on the unperturbed bath dynamics. Using the additional 
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assumption that the bath response is of the same kind as 
that studied in earlier publications, we derive the central 
theoretical result of this paper. In Section III we discuss 
the error associated with the basic assumption of Section 
II. In Section IV we use the central theoretical result of 
this paper to predict the resulting, non-canonical, equi- 
librium. In Section V we check numerically the effect of 
the error discussed in Section II. Finally, in Section VI 
we make a balance of the results obtained in this paper. 



II. THE PROJECTION METHOD 

In the case of no feedback, namely when A = 0, it is 
shown |[o| that an immediate benefit of the dichotomous 
choice is that the projection operator method allows 
us to express the dynamics of the variable of interest, x, 
in terms of an exact and simple diffusion-like equation of 
motion. When the feedback is included, unfortunately, 
the projection method does not produce a simple equa- 
tion of motion, and delicate assumptions must be made 
if we want to keep the elegance and simplicity of the ear- 
lier treatment. The main purpose of this section is to 
discuss these delicate assumptions. The adoption of the 
projection method yields for the part of interest Ppr 
of the total distribution the following equation: 



d_ 

at 



P PT (x, w, t) = J [P[~^] exp[QLQ(t - t')] ■ 

Q[-^ + n-A 2 x)]Pp T (x, t w, f )} dt'. (2) 



As usual, we assume the bath to have an equilibrium 
distribution satisfying the condition: 



rop eq (£,w) = 0. 



(3) 



This dictates the choice of the projection operator P: 

Pp T (x, £, w; t) = a(x, t)p eg (t;, w), (4) 

where o~{x, t) is obtained tracing the total distribution 
Pt{x, £, w, t) over the irrelevant degrees of freedom £ and 
w. Note that Eq.(|^) is an exact equation provided that 
the initial condition is given by: 



p T (x, £, w; 0) = a{x; 0)p eq ((, w). 
For simplicity we assume that: 

t(K)=t + Kt 1 . 



(5) 



(6) 



We carry out our calculations setting the condition: 
K = — A 2 x. This is equivalent to assuming a form of lin- 
ear response to external perturbation in agreement with 
Refs. [p2|-p^[ . We make now the assumption that, in spite 
of A 7^ 0, the exponential operator exp(Q LQ(t — t') ap- 
pearing in Eq. (0) only depends on the unperturbed op- 
erator Tq, a property that, as earlier pointed out, is valid 



only in the free diffusion case. We also use Eq.(||). Un- 
der all these approximations, we are allowed to rewrite 
Eq.(|) as 

d - 

—Pp T {x, f , w, t) = A(p T (x, f , w, tj) + B(p T (x, £, w, t)), 

(7) 

where 

A(p T (x,^w,t)) = — dt'|p[£]exp[ro(i-0]- 

Q[Z]PM*>Z>vJ)} ( 8 ) 

and 

B(pr(x^,w,t))=A 2 -^xJ^ dt'{P[£]eMh(t-t')] ■ 

<Hfi]PpT(x,£,W,i')}. 

(9) 

In conclusion to derive this result we have used a major 
assumption that will be referred to as assumption (^.This 
assumption can be expressed as follows: 

assumption (i). We assume that the exponential opera- 
tor appearing on the r.h.s of Eq.(Q) only depends on the 
unperturbed bath dynamics. 

In the case where no feedback process is considered 
|Tof there is no error associated with this assumption, 
since this is shown to be an exact consequence of the 
dichotomous nature of the variable £. It is not so in 
the more general case of this paper. We shall devote 
Section III and the numerical treatment of Section V to 
assessing the consequences of the error associated with 
this basic assumption. Adopting the formalism of the 
response theory jj] we rewrite Eq.(^) in the form: 



d d 2 
— <r(x,t) = < > eq J o dt'^(t-t')-Q^^,t' 



A 2 <£r 1 > eq J dt'C{t-t')—xcr(x,t'), (10) 



where 



and 



< £ 2 >eq 



C(t) 



< £ exp(f o t) f i > 



< <i 



>eq 



(11) 



(12) 



This result has been obtained by evaluating the diffusion 
term at the zero-th order in the feedback interaction, 
and considering, in agreement with the linear response 
criterion Q, the first non- vanishing contribution, pro- 
portional to the friction. It is worth remarking that the 
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correlation function of Eg. (|1 1|) affords the most conve- 
nient way of defining the microscopic time r mentioned 
in Section 1. This is given by: 



$t{t)dt. 



(13) 



We now have recourse to the second approximation on 
which our crucial theoretical results rests. This assump- 
tion will be referred to as assumption (ii) and can be 
expressed as follows: 

as sumption (ii). We assume that the function C(t) has 
a finite time scale. 

This assumption is dictated by the theoretical and nu- 
merical conclusion of the earlier work of Refs. [p^|-p^[. 
This assumption cannot be mislead as a property of or- 
dinary statistical mechanics. Actually, this assumption 
means a deviation from ordinary statistical mechanics, 
which, as shown in Ref. B, would imply 



C(t) = 



(14) 



In the case where the correlation function (t) is not in- 
tegrable and the correlation time of Eq. ( |l3| ) diverges, the 
condition of Eq.(p"4|) would imply a field of finite inten- 
sity to produce a current of infinite intensity fl2|| . The 
numerical calculations show that this striking physical 
condition is not realized fl2|| , thereby implying a viola- 
tion of Eq. (14). This violation, in turn, is due to the 



fact that the function C(t) has a finite time scale even 
when the function does not. Note that under the 

assumption that the function C(t) has a finite time scale 
it is possible to define 



7 = A < £Ti > £ 



dt'C{t'). 



(15) 



In conclusion, we obtain the following important equa- 
tion: 

d r* d 2 

— a(x,t) =< > eq J dt'<S>^t~t')—a{x,t r ) 

d 

+'Y-Q^xa(x,t). (16) 

This result rests on both assumption (i) and (ii). How- 
ever, it is evident that special attention must be devoted 
to the first assumption. In a sense the validity of assump- 
tion (ii) has already been assessed by the theoretical and 
numcr ical work of Refs. The validity of assump- 

tion (i), on the contrary, requires further discussion. This 
will be done in Sections III and V. 

Before ending this Section, we want to remark that 
in the special case where the condition of Eq.(p^) ap- 
plies, the important result of Eq. (|Io| ) becomes very simi- 
lar to the Fokker-Planck type of equation recently found 
by the authors of Ref. [[l5| . These authors pointed out 
that an equation of this kind shows that anomalous dif- 
fusion can be compatible with Boltzmann statistics. We 



note that this conclusion does not apply to the case of 
super-diffusion under study in this paper, because, as we 
have seen, Eq.([l4]) does not apply. In the subdiffusional 
case studied in Ref. M] , however, there are no compelling 
reasons leading to the breakdown of Eq. ( 14 ) thereby leav- 
ing open the possibility that the subdiffusional condition 
is compatible with Boltzmann statistics. 



III. NO INTERFERENCE BETWEEN FREE 
FLUCTUATION AND DISSIPATION:TIME 
EVOLUTION 

To a first sight, one might be led to think that Eq.(|l6|) 
is equivalent to the Langevin-like equation: 



j(t) = -7 *(*)+£(*). 



(17) 



supplemented, of course, by the set of equations neces- 
sary to determine the time evolution of the dichotomous 
variable In this Section we show that Eq.(|l6|) is 

not identical to the equation of motion for o~(x, t) gener- 
ated by Eq. (|l7|) . This will help us to estimate the error 
affecting the main prediction of this paper about the con- 
dition of equilibrium established by the feedback on the 
generator of fluctuation without time scale. 



A. Second moment time evolution 

In Section IV we shall point out that Eq. (|l7|) implies 
that throughout system's time evolution the trajectory 
x(t) departing from the initial condition a;(0) = never 
leaves the interval [—W/j,W/'y}. This property means 
that the second moment of the distribution is kept finite 
at all time and can never exceed the maximum value 
(W/^) 2 . Here we show that, on the contrary, the second 
moment of the distribution driven by Eq.(16) diverges for 
t — > oo. 

Using Eq.(ll^) we get: 



^<x 2 (t)>=jj dxx 2 -^-[xo-(x,t)} 

,--oo ft op, 

+ <e 2 >/ dxx 2 dt'^> i (t~t')—a(x,t / ). (18) 

Joo JO VX 

Using the method of integration by parts it is shown that 
Eq.pl) yields: 



d_ 

dt 



< x 2 (t) >= -27 < x 2 {t) > +2 < £ 2 > [ <$>(t')dt'. 

Jo 

(19) 



Note that the first term on the r.h.s. of Eq.(|l9|) can be 
derived from the first term on the r.h.s. of Eq.(^) via 
integration by parts, provided that the decay of the func- 
tion a(x,t) for |x| — > 00 is faster than l/|x| 3 . This means 
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that the distribution a(x,t) cannot be a Levy process at 
any finite time t > 0. We know that at 7 = the dif- 
fusing distribution is in fact a Levy process with ballistic 
peaks signalling the presence of a propagation front jio) , 
thereby ensuring the validity of the method of integra- 
tion by parts. It is plausible to assume that the action 
of a dissipation process makes the distribution spreading 
less intense, so favoring rather than opposing the method 
of integration by parts. 

The solution of Eq. ( |l9| ) is given by: 



< x 2 (t) >- 



<e> 
1 



$ ? (f - t')[l - exp(-2 7 <')]<&'• 



(20) 



Let us adopt for the correlation function the choice: 



(21) 



where T is the mean waiting time in a state of the ve- 
locity. In fact, as a consequence of the one-dimensional 
assumption we are allowed to use the relation JlCJ] : 



(22) 



where ip(t) is the distribution density of sojourn times. 
Plu ggin g the analytical form of Eq. |2l] into the r.h.s. of 
Eq. E0j and making a time asymptotic analysis, we get: 



and 



lim < x 2 (t) >« t 1 - 13 , 7 > 

t— »oo 



lim < x 2 (t) >« t 2 ~P, t » 1. 
7^0 



(23) 



(24) 



B. Exact equation of motion for a(x,t) 

We note that the use of the same projection method 
as that applied in Section II to the dynamic system de- 
scribed by Eq.(|l7|) yields 

§- t *(x,t) = 1 ^xa(x,t) + ^J*dt>{* d t-tr 

eM ^ x{t -t')]^-a{x,t')}. (25) 

We immediately see that the same approximation as that 
applied to Eq.(0), namely the approximation of neglect- 
ing the influence of the feedback on the memory ker- 
nel, makes Eq. ( p5| ) identical to Eq. (|l6|) . Consequently, 
the numerical treatment of Eq. ([17]) is expected to depart 
from the prediction of Eq. (|l^) and the amount of this de- 
parture can be used as a way to establish the error caused 



by assumption (i) in the derivation of Eq.jjlq), which is 
the central result of this paper. 

Eq.(p5f) can be used to derive an analytical expression 
for the second moment time evolution. The Taylor series 
expansion of the exponential operator on the r.h.s. of 
Eq.(^5|) and the use of integration by parts yield: 

d_ 

Of 



< x 2 {t) > +2 7 < x 2 {t) > = 2 < £ 2 > 



/ <& ? (i')exp(- 7 0< 
Jo 

(26) 

which, in turn, yields the following time evolution: 

< x 2 (t) >= 2 < £ 2 > cxp(— 2jt) / exp(27£')- 

Jo 

[ $^(r)exp(-7r)drdt / . (27) 
Jo 

It is worth remarking that the general expression for the 
asymptotic value of the second moment is: 



< x 2 (oo) >= 



<e> 



1 Jo 



exp(-7*')(2i' . (28) 



We see that the asymptotic value for the second mo- 
ment, as it must be, is finite and in the special case of 
Eq. ( pl| ) the analytical expression for the second moment 
at £ — 00 is: 



r(l-/3, 7 /?T) 



7 



2-/3 



(29) 



In conclusion, we see that assumption (i) produces the 
seemingly unacceptable effect of making the asymptotic 
second moment diverge, whereas the exact equation of 
motion yields a second moment which is always finite and 
gets at equilibrium the finite value predicted by Eq. (^9|) . 
The discussion of Sections IV and V will explain in which 
sense the error associated with assumption (i) does not 
invalidate our main conclusion that the final equilibrium 
distribution is of Levy kind. 



IV. EQUILIBRIUM PROPERTIES 

The fact that the second moment does not converge 
to a finite value is a consequence of the central approxi- 
mation yielding Eq. (|l6|) . This does not conflict with the 
possibility that for t — > 00 the distribution approaches 
asymptotically a time independent shape. Using the re- 
cent results of Refs ]icj and |ll| it is shown that the 
Fourier transform of Eq. ( |l6| ) obey the following time evo- 
lution equation: 

-a(k, t) = -b\k\ a a(k, t) - lk—a(k, t), (30) 



4 




< ar oo > 







FIG. 1. Equilibrium distribution Vix) = <r(a;,oo) as a 
function of x. The Levy distribution obtained with the 
anti-Fourier transform of Eq. ( |3l| ) (dashed line) fits very well 
the distribution middle. The system parameters used are, 
7 = 1CP 4 , T — 50, f3 = 0.6. The numerical conditions are: 



number of trajectories = 10 4 , observation time 
interval = 80 and b = 8.4 • 10 4 . 



2 ■ 10 J , bin 



where a = 1 + (3 and & is a positive constant (see Eq.(21) 
of Ref [^6|). In particular, the paper of Ref shows 
that in the asymptotic limit the contribution to free dif- 
fusion becomes indistinguishable from that originally in- 
troduced by Seshadri and West jl?]] . Equation ( J30| ) yields 
the following equilibrium distribution: 



&(k, oo) = exp f — I /cl"^ 

V a 7 / 



(31) 



which, in turns, according to p8| , coincides with the equi- 
librium distribution corresponding to the equation of mo- 
tion: 

d 



dt 



X(t) = -J X (t)+T){t), 



(32) 



where ?y(i) is an uncorrelated noise, with probability dis- 
tribution p(j]) , obeying Levy statistics, and thus defined 
in the Fourier space by: 



p(k) 



dr\ exp(— ikrj)p(j]) = exp(— b \k\^), (33) 



where < fi < 2. 

It must be pointed out that Eq.(|32|) does not coincide 
with Eq.(|l7|). In the case of very weak friction they do 
in the sense that will be illustrated in Section V. 



V. NUMERICAL RESULTS 

The numerical results of this Section are based on the 
numerical treatment of Eq. ( |l7| ) , and consequently on the 
numerical implementation of 



FIG. 2. Comparison between computed and theoretical dis- 
tribution variance. The theoretical prediction used is Eq. 

The crosses indicate the numerical result and the lines 
the theoretical prediction. Each curve has been obtained 
keeping 7 constant. From the bottom to the top curve the 
values of 7 are: 10" 2 , 5 ■ 10" 3 , 10" 4 , 10" 5 . 



"(*) = 



exp[ 



-lit - t')W) dt' + x(0) exp(-Tt) 



(34) 



The fact that the variable £ is dichotomous with the cor- 
relation function of Eq.(|2l|) naturally leads us to adopt 
the same numerical approach as that used in Refs. p0| , ^9| . 
This means that two random number generators are used. 
The first results in random number homogeneously dis- 
tributed in the interval [0, 1]. With a proper non-linear 
deformation this is made equivalent to a random gener- 
ation of waiting times with the distribution of Eq. (pi]) . 
This is the way we adopt to build numerically the time 
evolution of x(t). We also set an initial condition fit- 
ting the crucial condition of Eq.(|^) and we make the 
trajectory run for times larger than 2O/7. We evaluate 
the same way 10 4 trajectories, then these trajectories are 
recorded in a bin. 

In Fig 1 we see a sample of the resulting equilibrium 
distribution with (3 = 0.6, T = 50, W = 1, 7 = 10~ 4 , 
spanning from x — —W/j to x = W/"/. Fig. 1 is a crystal 
clear illustration of what we mean by statistics of Levy 
kind. We see that the equilibrium distribution is trun- 
cated and that at the borders two sharp peaks emerge. 
These sharp peaks are an equilibrium reflection of the 
peaks already revealed by the numerical treatment of free 
diffusion [jlO 20 1. However, the distribution enclosed by 
these peaks is shown to fit very well the Levy distribution 
predicted by the theory of Section IV. 

Fig. 2 is devoted to the comparison between the the- 
oretical prediction of Eq.(j2^) and the result of our nu- 
merical treatment. The agreement between theory and 
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numerical calculation is extremely good, not only in the 
Gaussian regime j3 > 1, but also in the regime f3 < 1 
up to (3 ~ 0.5 . Significant discrepancies between theory 
and numerical results can be found in the region close 
to (3 — 0, probably as a consequence of the fact that 
with a finite number of trajectories the peaks, which are 
expected to give significant contribution to the second 
moment, cannot be satisfactorily reproduced. 

In conclusion, these numerical results prove that we are 
in a full control of the error caused by approximation (i). 
The markedly Levy character of the equilibrium distri- 
bution in the sense here illustrated with the help of Fig. 
1 supports our theoretical prediction: The main feedback 
effect is the emergence of the non-canonical statistics of 
Ref. & 



VI. CONCLUDING REMARKS 

Which is then the interest of our results? We think 
that the interest of them lies on this: This paper forces 
us to change the conventional perspective concerning the 
microscopic foundation of the canonical statistical be- 
havior. Some years ago, the findings of Zhu and Robin- 
son |2l|] have been criticized by Keirstad and Wilson p^ ] 
with arguments which are a nice example of the con- 
ventional wisdom. Let us see why. Zhu and Robinson 
plf had detected significant deviations from the canoni- 
cal Maxwell velocity distribution, in a physical condition 
characterized by a system of interest very fast compared 
to its thermal bath. This condition seems to be related 
to that considered in this paper where the dynamical sys- 
tem playing the role of bath is in fact so slow as to break 
the condition itself of time scale separation. The reaction 
of the scientific community, of which the authors of p^ j 
are a significant example, has been that the non canoni- 
cal behavior detected numerically by Zhu and Robinson 
pl[ is an artefact of numerical inaccuracy and limited 
computation time. This paper shows, on the contrary, 
that the opposite condition might apply, namely, that 
ordinary rather that anomalous statistics might be the 
result of numerical inaccuracy. We know that the round- 
off errors are equivalent to the influence of fluctuations 
of a given intensity e. The larger the computer accuracy, 
the smaller is the intensity of the equivalent fluctuations. 
On the other hand, we know [^3| that the effect of these 
fluctuations is that of changing the correlation function 
of Eq.(ll) into the correlation function <E>|(i) related to 
the original by 



(35) 



with tc proportional to e 5 and 5 being a positive coef- 
ficient, of the order of unity, determined by the micro- 
scopic dynamics under study p3[ . It is evident that at 
times t > tc the Markov approximation is valid, and as 
an effect of it the non standard equation of ( [Toj ) becomes 
identical to a conventional Fokker-Planck equation. The 



non conventional equilibrium of ( j25| ) is a time asymptotic 
property, and at any given time t > > I/7 we can produce 
a transition from the regime of non-ordinary statistics to 
a regime of canonical Gaussian equilibrium by increas- 
ing the intensity of the parameter e, so as to realize the 
condition I/7 > Tc- 

Finally, we want to stress a problem worth of future in- 
vestigation. This has to do with the increasing attention 
devoted to the non-extensive thermodynamics of Tsallis 
pi] p6[ . Non-extensive thermodynamics means that the 
deviation from the canonical equilibrium distribution is 
not more perceived as a violation of statistical mechan- 
ics. This is a very valuable aspect of this research work 
pi] p6| . In fact, as a result of the interest that Tsallis's 
non-extensive statistical mechanics is raising, a deviation 
from the ordinary prescription, of the kind earlier men- 
tioned, would be judged these days as a possible manifes- 
tation of non-extensive thermodynamics triggered by the 
long-range correlations of the dynamical system under 
study, rather than a consequence of numerical inaccu- 
racy. 

However, the arguments of this paper show that under 
the specific form here adopted to establish a fluctuation- 
dissipation process in a case of dynamics without time 
scale, the basin of attraction for equilibrium distribu- 
tion is given by Levy statistics. It is interesting to point 
out that Levy statistics share with Tsallis statistics the 
power law behavior of the distribution tails. However, 
the central part of the Levy distribution significantly de- 
part from the generalized canonical distribution of Tsal- 
lis. In an earlier paper |27| it has been shown that the 
adoption of Tsallis' non-extensive thermodynamics natu- 
rally leads, via entropy maximisation under a proper con- 
straint, to a transition probability with an inverse power 
law decay at large distances. By repeated application 
of this kind of transition, as a consequence of the Levy- 
Gnedenko theorem |^8| the diffusion process is attracted 
by the basin of Levy statistics. In the case of extremely 
weak friction, equilibrium is reached as a result of a very 
large number of elementary transitions, and this is prob- 
ably the main reason why eventually the resulting statis- 
tics is of Levy rather than Tsallis kind. 
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